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Abstract 

We compute the response and the angular pattern function of an interferometer for a 
scalar component of gravitational radiation in Brans-Dicke theory. We examine the prob- 
lem of detecting a stochastic background of scalar GWs and compute the scalar overlap 
reduction function in the correlation between an interferometer and the monopole mode 
of a resonant sphere. While the correlation between two interferometers is maximized 
taking them as close as possible, the interferometer-sphere correlation is maximized at a 
finite value of / x d, where / is the resonance frequency of the sphere and d the distance 
between the detectors. This defines an optimal resonance frequency of the sphere as a 
function of the distance. For the correlation between the Virgo interferometer located 
near Pisa and a sphere located in Frascati, near Rome, we find an optimal resonance 
frequency / ~ 590 Hz. We also briefly discuss the difficulties in applying this analysis to 
the dilaton and moduli fields predicted by string theory. 



1 Introduction 



A number of interferometers for gravitational wave (GW) detection are presently under con- 
structions and are expected to be operating in the next few years. In particular, VIRGO is 
being built near Pisa, the two LIGO interferometers are being built in the US, GEO600 near 
Hannover, and TAMA300 in Japan. These interferometers are in principle sensitive also to a 
hypothetical scalar component of gravitational radiation. Scalar GWs appear already in the 
simplest generalization of General Relativity, namely Brans-Dicke theory, whose action reads 



Sbd = / f/! - r v 



ipR - ^v'VV uV 



+ 5 ma tter > (1) 



with ip the Brans-Dicke scalar. The coupling of matter with gravity, 5 ma tter; is dictated by 
the equivalence principle. In order to avoid conflict with solar system experiments, one must 
take |wbd| greater than approximately 600. 

At a more fundamental level, various scalar fields with interactions of gravitational strength 
come from string theory. A universal example is the dilaton $; the low energy effective action 
of string theory, in the graviton-dilaton sector, reduces to the first term on the right-hand 
side of eq. ([!]), with the identifications ip = (16ir/a')e~ 2 ® (where a' is the string tension) 
and lubd = — 1- To avoid conflict with experiments, it is expected that the dilaton will get 
a mass from non-perturbative mechanisms Q, or that it decouples from matter with the 
cosmological mechanism proposed in [^]. Furthermore, various scalar fields (moduli) appear 
when compactifying string theory from ten to four dimensions. Their number and couplings 
are strongly dependent on the specific compactification used. 

In this paper we investigate whether it is possible to search for such scalar particles using 
the GW interferometers under construction, as well as the resonant spheres which are under 
study. We start from the Brans-Dicke theory and, in sect. 2, we discuss the response of an 
interferometer to a GW with a scalar component: in particular, we find that such a scalar 
component creates a transverse (with respect to the direction of propagation of the GW) 
stress in the detector; we compute the phase shift Aip measured in the interferometer and 
derive the angular pattern function, i.e. the dependence of the signal on the direction (6, </>) 
of the impinging GW (see figure |l|). We find A(p oc sin 2 9 cos 2<fi. We also show the physical 
(and formal) equivalence of two different gauges used to describe scalar radiation. 

In sect. 3 we consider the detection of a stochastic background of scalar GWs. In this 
case it is necessary to correlate two different detectors. We give a general treatment of 
the computation of the overlap reduction function T(/) for generic detectors in the scalar 
case: such a function represents a measure of the correlation between the signals of the two 
detectors and depends on the frequency of observation / and on the type of detectors one 
uses, as well as on their location and relative orientation. Similarly to what has been done in 
ref. [14] in the case of the + and x components, we "factorize out" from r(/) the response 
tensors Dij of the detectors, which summarize the whole information about the type of the 
detectors and their orientation in space; next we compute explicitly the remaining part of 
r(/), that is the dependence on the frequency and the location of the detectors. The result 
is thus completely general and it is applicable to any given pair of detectors. 

We then examine in particular the correlation between the VIRGO interferometer and a 
resonant sphere, as the prototype which is presently under study in Frascati, near Rome || ||, 
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||, ||, (7|. Similar correlations are also in principle possible between LIGO and the TIGA reso- 
nant sphere located in Louisiana || . We compute the interferometer-sphere overlap reduction 
function and we find that, contrarily to what happens in the correlation of two interferome- 
ters, the correlation is not optimized when the detectors are as close as possible (compatibly 
with the constraint of decorrelating local noises), but instead there is an optimum non zero 
value of the product of the distance between the detectors and the resonance frequency of 
the sphere. For the distance between the VIRGO interferometer in Cascina, near Pisa, and 
a resonant sphere in Frascati, near Rome, we find an optimum correlation if the resonance 
frequency of the sphere is / ~ 590 Hz. This value of / is quite interesting because is in 
the range where interferometers achieve their highest sensitivities and at the same time is 
comparable to realistic values for the resonance frequency of the spheres which are presently 
under study. Actually, the TIGA prototype has its first resonant mode at 3.2 kHz, but hollow 
spheres, which are presently at the stage of preliminary feasibility studies ||, depending on 
the material used and other parameters, could have a resonance frequency between 200 Hz 
and 1-2 kHz, with a bandwidth of order 20 Hz. 

In sect. 4 we discuss how the interaction with the detectors is modified by a very small mass 
term for the scalar field. This is partly motivated by the desire to examine the perspective 
for detection of the string dilaton and moduli. We find that, in presence of such a mass term, 
the stress induced in the detector by the scalar wave is not anymore purely transverse, but 
has also a longitudinal component of relative amplitude m 2 /w 2 , if m is the mass of the scalar 
and uj is the frequency of the GW. 

In sect. 5 we will then briefly discuss some difficulties in applying our analysis to these fields 
predicted by string theory. A number of technical details are collected in the appendixes. 

2 The response of the interferometer to scalar gravitational 

waves 

2.1 Computation in the gauge eg = diag(0, 1,1,0) 

In this section we compute the phase shift measured in the interferometer when a scalar GW 
is coming from an arbitrary direction. There are of course different possible gauge choices 
(i.e. coordinate transformations) for representing plane wave solutions of the equations of 
motion of Brans-Dicke theory with both spin 2 and spin components. We first consider a 
gauge choice that, for a wave propagating in the +z direction, brings the metric perturbation 
in the form 



V(* ~z) = AM(jt - z) e<+> + A^\t - z) e M + *(i - z) eg 



(2) 



where e^' x,s 



are the polarization tensors 
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The choice of gauge that brings the plane wave solutions of the equation of motion of Brans- 
Dicke theory into this form is discussed in ref. [[h], ||] and, for comparison with later results, 
we recall the main points of the derivation in app. A. Under rotations around the z axis, 
it is straightforward to verify that e ^ db i have helicities ±2, while has helicity 0. 
Thus the term A^^ eju) + e)£) describes ordinary gravitational waves with + and x 

(s) 

polarizations in the transverse traceless gauge, and the term $ e^l, describe a scalar GW, 
characteristic of the theory that we are considering. 

To compute the response of the interferometer we start with the geodesic equation for a 
free-falling mass, 

sr + x a = , (5) 

where 

r «/3 = (dahpv + dph av - d u h a/3 ) (6) 

are the linearized Christoffel symbols and ( ) = — denotes derivation with respect to proper- 
time of the mass. In the gauge (|J) we have Tq = and therefore, if the mass is initially 
at rest, it will remain at rest also when the wave arrives, and the proper time r is the same 
as the time variable t. In other words, choosing the gauge (Q) means that we automatically 
choose the reference frame whose coordinates are comoving with free-falling masses. 

The situation is then perfectly analogous to the ordinary gravitational waves in the trans- 
verse traceless gauge: the coordinates of the test masses are unaffected by the gravitational 
wave, but physical ( i.e. proper) distances between masses are influenced by the wave. In our 
case the line element is: 

ds 2 = - dt 2 + (Sij + hij(t - z)) dx { dx j = 

= - dt 2 + dz 2 + (1 + $(t - z)) (dx 2 + dy 2 ) , (7) 

where we have restricted ourselves to a purely scalar wave and we use units c = 1; as usual 
we suppose that the wavelength of the scalar gravitational wave is much larger than the 
distance between the test masses, which for the interferometer are the two mirrors and the 
beam-splitter; we take the latter at the origin O of the coordinate system (note again that in 
this gauge if the beam-splitter is at rest in the origin before the waves comes, it will always 
remain there; this will not be true in the gauge considered in the next subsection) and so 
its frequency / is much smaller than 1/To, where To is the time a laser-beam takes to travel 
from O to the mass; under this assumption, we can consider the amplitude — z) frozen at 
a value <3?o- So, the proper distance of the mass of coordinates X, Y, Z from the origin O is: 

L = v /(l + <D )(X2+y2 ) + Z 2 . (8) 

The physical interpretation of @ is clear: the wave acts only on the coordinates that are 
transverse with respect to its direction of propagation. If we denote by L = VX 2 + Y 2 + Z 2 = 
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To the proper distance that the mass has from the origin O before the wave arrives, then 
from eq. (H|) we see that, at first order in $o, 

ifZ = L=(l + ^ )L (g) 

if X = Y = => L = L . 

So, just like the ordinary gravitational waves with + and x polarizations, the scalar wave is 
transverse, not only in its mathematical description, eq. (||), but also in its physical effect. 

It is now easy to compute the phase-shift produced in an interferometer. We take an 
interferometer whose arms are aligned along the x and z axes and we consider first a scalar 
wave traveling in the +z direction. From eq. (^) we see that the time that the laser-beam 
takes in order to make N round-trips from the beam-splitter to the mirror in the x direction 
is 

T x = (l + i$ )iV2Lo, (10) 

while, for the beam traveling in the z direction, the travel time is unaffected by the scalar 
wave, 

T Z = N2L . (11) 

The phase-shift Aip is obtained (see e.g. [l^]) by multiplying the time difference (T x — T z ) 
by the angular frequency f2i aser of the laser; this leads to the result 

Aip = i$ • ip arm , (12) 

where (p arm = r2i aser iV2Lo is the phase that the laser-beam accumulates in TV round-trips. 

To compute the response of the interferometer for an arbitrary direction of propagation of 
the wave we recall that it is possible to associate to a detector a 'response tensor' D^ such 
that the signal (the phase-shift, if the detector is an interferometer) induced in the detector 



by a gravitational wave of polarization is proportional to Tr{D e} = eji [13, 14]; for 



an interferometer whose arms are along the u and v directions the detector tensor is 

= UiUj - ViVj . (13) 

The polarization tensor for a purely scalar wave traveling along a generic direction n is (see 
©) 

e\f(n) = 5ij - hihj , (14) 

and then the pattern function F s (h) = , which describes the interaction between a 

scalar wave propagating along n and an interferometric detector with arms along u and v, is 

F s (n) = -sin 2 6» cos 20 , (15) 

where 9 and <f> are, respectively, the polar and azimuthal angles of the versor n in the reference 
frame in which the x and y axes are co-aligned with u and v, (see figure [l]). Of course, if 
9 = 0, there is no phase shift because the proper length of the two arms is modified in the 
same way, and it cancels taking the difference. However the proper time that the light takes 
to make a round trip in each of the two arms separately is modified by the wave, and it is in 
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Figure 1: The definitions of the versors and angles discussed in the text. 



principle measurable. [] For a generic angle of incidence the cancellation does not take place 
and there is a phase shift. 

Combining (|l5| ) with ( |l2|) we get the phase shift for arbitrary direction of incidence of the 
scalar wave Jl6[ , 

A<p(0, (f>) = -~$0 • ^arm ■ Sm 2 9 COS 2(f) . (16) 

The angular sensitivity of the interferometer to scalar GWs is shown in fig. (||) , together with 
the standard angular pattern for the + and x polarizations. 



1 VIRGO takes its output from a set of 5 photodiods which give 5 optical lengths including the difference 
between the two arms, of course, and also the common mode, i.e. the sum of the two arm lenghts. However 
for the common mode the fluctuations in the laser power do not cancel, and the sensibility is quite limited; 
it is basically the same sensibility that could be obtained measuring the common mode with two separate 
interferometers, one made with one arms and the prestabilization cavity and the other with the other arm 
and again the prestabilization cavity. 
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Figure 2: The angular sensitivity of the interferometer for scalar waves (top), + polarization 
(bottom right) and x polarization (bottom left). The solid lines in the middle of the figure 
indicate the orientation of the arms of the interferometer. 
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2.2 



Computation in the gauge 



fJtU 



It is instructive to repeat the computation in a different gauge. Indeed, with a different 
gauge choice, we can bring the plane wave solution of the equations of motion into the form 



of eq. (Q) with effi and efi' unchanged, while 



so that, for a purely scalar wave, 



g^ = (1 + 



(17) 



(18) 



In app. B we prove this assertion and we find the coordinate transformation that allows to 
move from the previous gauge to this one.^ 

As in section 2.1, we consider a purely scalar gravitational wave traveling in the +z direc- 
tion and impinging on an interferometer whose arms are aligned along the x and z axes, and 
have a length Lq when there is no gravitational wave. We will denote quantities relative to 
the arms by subscripts X and Z. Again, we deal only with the case in which the frequency 
/ of the wave is much smaller then 1/Tq = 1/Lo, as discussed in the previous section. 

Consider the metric 



ds 2 



[1 + $)(-di 2 + dx 2 + dy 2 + dz 2 



(19) 



The laser-beam follows null geodesies (ds 2 = 0): so, for the beam that propagates along the 
x axis 

<lrLs,r ±1 ^x laser (t) = const ±t, (20) 



and for the beam along z 



dt 



dz\ a 



dt 



±1 =4> Z lager (t) = const ± t 



(21) 



The geodesic equation of motion for the mirrors and beam-splitter, in this gauge, is ]T 



where 



x(t) 

y(t) 

z{t) 

T{t) 

Kt- 



x 

JJO 

z + - z{t)) 
t + z{t), 



(22) 



z = 



t-z 



§{u)du . 



(23) 



2 This gauge has been used in ref. Jl^]. However, the authors did not realize that in this reference frame 
the beam splitter is not left at the origin by the passage of the GW, see below, and furthermore computed a 
coordinate-time interval rather than a proper-time interval, reaching the incorrect conclusion that the scalar 
GW has a longitudinal effect, and no transverse effect. Their expression for the phase shift and angular pattern 
then differs from ours. In this subsection we will confirm the results of section 2.1 , with a proper treatement 
of the scalar GW in this gauge. 
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(note that the solution for z(t) is an implicit equation). Therefore, in this gauge, the trans- 
verse coordinates are not affected by the passage of the wave, while the longitudinal coordinate 
is affected. Of course, this does not mean that the physical effect of the wave is longitudinal. 
The physical effect is determined looking at diffeomorfism invariant quantities, i.e. proper 
distances and proper times. 

Since in this gauge the x-coordinates of the mirrors and the beam-splitter are unaffected 



by the passage of the wave, from eq. (20) we see that the interval, in coordinate time t, that 



the laser takes for one round-trip in the X arm is 

T x = 2L . (24) 

Tx measures the time that the laser-beam takes in terms of the variable t: the beam leaves 
the beam-splitter at t = and comes back at t = Tx- However, this quantity is by definition 
not invariant under coordinate transformations, and in order to deal with interference at 
the beam-splitter, we must work in terms of the beam-splitter proper time; it is this latter 
variable that measures the physical length of the arms. Thus we call r(t) and zbs(^) the- 
proper time and z-coordinate of the beam-splitter at (coordinate) time t (with initial condition 
zbs(— 00 ) = 0). From eqs. (^), 

ZBs(t) = \ / $(u)du (25) 

J-OO 

, rt-z BS (t) 

r(t) = t + | / ${u)du. (26) 

J — oo 

The proper time interval that the beam takes to make a round-trip in the X arm is then 

tx = t(T x ) - r(0) = T x + ^ ${u)du ~ 

•^bs(o) 

~ T x + ~$a[Tx + zbs(0) - zbs(Tx)} ^ 

~ 2L (l + |l>o), (27) 

where we have used eq. (|24]), we have considered $ "frozen" at $0 an d we have kept only 
the first order in $ - Note that z B s(0) - Zbs(Tx) = O(&o). 

We now compute Tz and tz- Suppose that the beam leaves the beam-splitter at t = and 
reaches the mirror along the z axis at t = T\. Then 

2tos(0)+Ti = Z(ri), (28) 

where Z(t) denotes the z-coordinate of the Z mirror at time t (with Z{— 00) = Lq). Similarly, 
when coming back from the mirror the beam reaches again the beam-splitter at t = Tz = 
T\ + Ti such that 

Z(T 1 )-T 2 = z BS (T z ) . (29) 



Subtracting eq. (E9) from eq.pq) one has 

T Z = T 1 +T 2 = [Z(2i) - zbs(0)] + [Z{T X ) - z BS (T z )} . (30) 
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The equations of motion ( f22[ ) for the Z mirror and the beam-splitter read 

rt-Z(t) 



Z(t) 



<£(u)du 



oo 

*-^Bs(*) 



<&(u)du , 



which, substituted in eq. fl30|), give 

rT t -Z{T t ) 



T z = 2L + j 



^BS(O) 



Ti-Z(Ti) 
T z -z BS (T z ) 



<&(u)du 



(31) 
(32) 

(33) 



Because of eq. fl28|) , the first integral in eq. ( |33| ) is zero. The second integral, in the approxi- 
mation of a frozen wave and to first order in gives 



Ti-Z(Ti) 
Tz-z BS (Tz) 



$(u)du 



i$o [Ti - Z(Ti) -T z + z BS (T z )] ^ 
i$o[^o - - 2L + 0] = -|* 2L 



Eq. (|33|) thus becomes 



Following the same steps leading to eq. fl27|), we then obtain 

tz = t(T z ) - r(0) ~ T z (l + |*q) ~ 2L • 



(34) 



(35) 



Eqs. ( |27| ) and (p5[ ) agree with the result of sect. 2, see eqs. ( |10|) and (|il|), as they should. 
It is also apparent that the computation in the gauge used in the previous section is much 
simpler, due to the fact that in this case the coordinates of test masses are not affected by 
the passage of the wave, and only proper distances change. 



3 Stochastic backgrounds of scalar GWs 
3.1 General definitions 

The standard procedure described in [[lj, [l?], for detection of stochastic backgrounds of 
ordinary gravitational waves can be applied with minor modifications to the case of scalar 
GWs. A stochastic, gaussian and isotropic background of scalar GWs can be characterized 
by the quantity 

1 dp v 
Pc dlog/ 

where p v is the energy density associated to the background, / is the frequency and p c is the 
critical energy density for closing the Universe, 
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(Hq = 100 ho (km/sec) (1/Mpc) is the Hubble constant). The intensity of the background is 
expected to be well below the noise level of a detector: the detection strategy then consists 
in correlating the outputs of two (or more) detectors, located far enough so that local noises 
as the seismic noise are uncorrelated. One defines the quantity 



S 



-T/2 
-T/2 



dt 



-T/2 
-T/2 



dt' Sl (t) s 2 (t')Q(t-t') 



(38) 



where s a denotes the output of the a-th detector, T the total observation time and Q is a 
real filter function, that, for any given form of the signal, can be determined exactly in order 
to maximize the signal-to-noise ratio (SNR). We write s a (t) = h a (t) +n a (t), where h a is the 
signal induced by the background of scalar GWs (for an interferometer it is the phase shift 
computed in the previous section), and n a is the intrinsic noise of the a-th detector. Under 
the assumption of uncorrelated noises, the ensemble average (denoted by (•••)) of the Fourier 
components of the noise satisfies 



«(f) Mf')) = s(f-f')s ab stH\f\) 



(39) 



The functions Sn\\f\) are known as square spectral noise densities. The signal-to-noise ratio 
is defined as 



1/4 



SNRs [«^7s>F>J ' (40> 

where we have used the exponent 1/4 (instead of 1/2) to take into account the fact that S is 
quadratic in the signals (see (|38"|)). Optimal filtering (see |Tj], [l8|]) gives 



SNR 



2T 



1 sm U) 



1/4 



(41) 



The various quantities that appear in the SNR are defined as follows. 



The function Sh depends uniquely on the spectrum of the background, and not on the 
features of the detectors, and is related to by 



Shif) 



mi 



i 



(42) 



4vr 2 / 3 2 + 0->bd 

Note that in the case of ordinary gravitational waves the factor 1/(2 + u>bd) is absent 
The derivation of this result is given in app. C. 



• The function S n is defined as S n (f) = \J Sn\f) Sn'if) and therefore depends uniquely 
on the intrinsic noises of the detectors. 

• The function T(/) gives a measure of the correlation between the detectors, and depends 
on their relative position and orientation; it is defined as 

•dftfi 



:(2), 



r(/) 



47T 



(43) 



where d ■ m = x'" — x' > is the vector joining the detectors (d is their distance, rh a 
versor) and the (n) are the pattern functions of the detectors for scalar waves. 
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In the case of the correlation between two interferometers, it is conventional to define the 
overlap reduction function 7 by [19, 14 1 



T(f) 

7(/) = -f 2 ■ (44) 

where, for ordinary GWs 

FU=[^ E ^V^f (^aligned, (45) 

and the subscript means that we must compute F\2 taking the two interferometers to be 
perfectly aligned, rather than with their actual orientation. 

This normalization is useful in the case of two interferometers, since in this case, for 
ordinary GWs F\2 = 2/5, (while for scalar GWs F\2 = 4/15) and it takes into account 
the reduction in sensitivity due to the angular pattern, already present in the case of one 
interferometer, and therefore 7(/) separately takes into account the effect of the separation 
Ax between the interferometers, and of their relative orientation. With this definition, 7(/) = 
1 if the separation Ax = and if the detectors are perfectly aligned. 

This normalization is instead impossible when one considers the correlation between an 
interferometer and the monopole mode of a resonant sphere, since in this case F\2 = 0, as 
we will see below. Then one simply uses T(/), which is the quantity that enters directly 
eq. (fil]). Furthermore, the use of T(f) is more convenient when we want to write equations 
that hold independently of what detectors (interferometers, bars, or spheres) are used in the 
correlation. In the following we will always refer to T(f) as the overlap reduction function. 



3.2 The overlap reduction function for scalar GWs 

In this section we compute analytically the overlap reduction function of two generic detectors 
in the case of a background of scalar waves, generalizing the result for ordinary gravitational 
waves of ref. [Q. 

As one can see from eq. d43|), the overlap reduction function is obtained by averaging, over 
the possible directions of propagation of the gravitational wave, the product of the pattern 
functions of the two detectors, weighted with a phase that depends on the delay in the 
propagation from one detector to the other. The pattern function and detector tensors are 
related by 

F s (n) = A.eg^n) (46) 

As in the case of ordinary GWs, the response tensors D^' are normalized in such a way that 
the pattern functions Fs a \n) take 1 as maximum value, varying n. Then we can rewrite 



eq. (42) in the form 



r(/) = Dff r m {a,A) D® , (47) 



where a = 2irfd and 

iWa, A) = f^ 4 s) (n) e jJ(fi) e iM . (48) 
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As we discussed in the previous section, the polarization tensor for a scalar wave traveling in 
the direction n is 



Xj,Xj + ViVj — Oij 



(49) 



where x and y are versors perpendicular to n and to each-other (see fig. [T]). The difference 
between the scalar case and the traditional one resides in T^m, because in the latter case the 

polarization tensors e^~' x ^(n) take the form: 



elf (&) 



(50) 



We now proceed to compute the tensor r^jy. Note that the result we find is absolutely general, 
independently on the type of the detectors used in the correlation, since the informations on 
the detectors is summarized in the response tensors D^' and not in T 



ijkl ■ 



Since the tensor JT^-jy is symmetric in i *-* j, in k I, and in (ij) <-> (kl) (see eq. (|48|)), it 
can be written in the most general form as 



where 



T ijk i(a,m) = A(a)A ijk i + B(a)B ijH + C(a)C ijM + D(a)V ijkl + E(a)£ ijH 



Bijki = SikSji + SuSjk 

Cijki = Sijrhkrhi + 5 k irrnrhj 

Vijkl = mimjm k mi 

£ijki = 8 ik rhjrhi + 8j k rhirhi + 5jirhirh k + 5urhjm k 



(51) 



(52) 



(the tensor ei km m m eji n rh n + €ii m m m ej kn rh n is a linear combination of these). 

To determine the coefficients A(a), B(a), E(a), one has to contract Tijki (written 



both in the form (48) and in the form (pl[)) with the five tensors (|52| ) and to calculate the 
corresponding scalar integrals, which can be done in terms of spherical Bessel functions; one 
thus obtain a linear system with five equations, whose solution is 



" A ' 




B 




C 
D 


. E . 





a 2 • jo (") - 2a • j x (a) + h ( a ) 

.72(a) 

-a 2 • jo (a) + 4a • ji(a) - 5 • j 2 (a) 



a 



3o(a) 



10a • ji(a) + 35 • .72(a) 
a • ji(a) - 5 -.72(a) 



(53) 



sin a 



a 



where the jg are the spherical Bessel functions 

Jo («) 
Ji(a) 
.72 (a) 



io(a) — cos a 
a 

o Ji(a) 



a 



Jo (a) 
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Figure 3: Relative position of the sphere and interferometer. 



The overlap reduction function for any specific two-detectors correlation can now be obtained 
using this general form for r^jy and then plugging into eq. (^7|) the appropriate detector 
tensors. 

A particularly interesting correlation is the one between the VIRGO interferometer and 
a resonant sphere. A resonant sphere has 6 different detection channels: 5 corresponding 
to the quadrupole modes with i = 2, m = 0, ±1, ±2, and 1 corresponding to the monopole 
mode with I = 0, m = 0. The monopole mode is especially interesting because it cannot be 
excited by ordinary spin 2 GWs. The response of a resonant sphere to scalar GWs has been 
computed in ref. ||. 

To correlate an interferometer, with arms along u and v, with the monopole mode of a 
resonant sphere, we use the response tensors 

£)W = u<g>u-v(g)v (54) 

= Ii 3 , (55) 

and we get 

r(/) = (sin 2 #cos2^) j 2 (2irfd) , (56) 

where 6 and cj> are the angular coordinates of the resonant sphere (i.e. of m) with respect to 
the "natural" reference frame of the interferometer (see figure ||) . 

Figure |] shows the frequency dependence of this overlap reduction function for a particular 
choice of the distance d corresponding to the distance between the VIRGO site and Frascati. 
The most striking aspect, compared to the overlap reduction functions of interferometer- 
interferometer correlations, is that the correlation vanishes if the distance between the de- 
tectors, d, goes to zero. This is the opposite of what happens for two interferometers, where 
instead the correlation is maximum for coincident detectors, so that the ideal strategy in 
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Figure 4: The behavior of j2(27r/d) against the frequency / (Hz) for d = 270 km, which is 
approximately the distance between the location of VIRGO and Frascati (resonant spheres). 



that case would be to have two interferometers as close as possible, compatibly with the 
fact that local noises, like the seismic noise and local electromagnetic disturbances should be 
uncorrelated (it is believed that this implies a minimum distance of at least a few tens of 
kms). In this case, instead, a sphere coincident with an interferometer makes for a totally 
ineffective correlation. This is easily understood, because the monopole mode of the sphere 
has, obviously, a constant angular pattern, while an interferometer has a dependence ~ cos2(/> 
which over the solid angle integrates to zero. Parentetically, this shows that in this case the 
function 7(/) cannot be defined, since it should be defined dividing T(/) for the value for 
coincident detectors, which is now zero. However, T(f) is the quantity that enters the SNR. 

Of course, as d — > oo the overlap reduction function again vanishes, and therefore there is 
an optimal, finite value for the product of the distance d and the frequency /, determined by 
the maximum of j2(2nfd). Approximately, this maximum is reached when 

f ^ d ^ ~ 1 . (57) 



590H Z y V 270km 

We have taken the reference value 270 km, which corresponds to the distance between the 
locations of VIRGO (Lat. N 43.63, Lon. E 10.50 deg) and Frascati (Lat. N 41.80, Lon. E 
12.67 deg). Since the resonant sphere is a narrow band detector, compared to the interfer- 
ometer, this means that the optimal correlation is obtained with a sphere that resonate at 
approximately 590 Hz. We see from fig. (|j) that the sensitivity of the overall correlation is 
very strongly affected by the value of /, and a value of / ~ 1 kHz, while keeping the sphere 
in Frascati, can easily result in the loss of more than one order of magnitude in the SNR 
defined in eq. (41). Note also that the minimum detectable value of h^O,^ is quadratic with 
the SNR and therefore an inappropriate value of the resonant frquency can easily result 
in loosing two or three order of magnitude in h^Q^. 

Conversely, if for technological reasons the frequency of the resonant sphere has to be fixed 



at a different value, eq. (|57 ) can be used to obtain the optimal distance between the two 
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detectors. 



4 The interaction of the interferometer with a very light scalar 

Till now we have considered only massless scalars; it is interesting, having in mind the 
extension to situations motivated by string theory, to consider the effect of a small mass term 
on the response of a detector. 

Let us first understand what it means 'small' in this context. We still want to treat the 
scalar as a classical wave that acts coherently on the detector. This implies first of all that 
m <C where L is the characteristic size of the detector. For a sphere of size of order 1 
meter, this gives approximately m < 10~ 6 eV, while for a km sized interferometer m < 
eV. In both cases, however, there is a stronger limitation that comes from the fact that, if 
we want to detect these scalars with a detector that works at a frequency fo = 1 kHz, or 
ujq = 2ir kHz, we must require that the frequency of the massive scalar, ui m = y/k 2 + m 2 , is 
of order of loq. Since of course u m > m, we get 

m < 2yrkHz ~ 4 x l(T 12 eV. (58) 

For these ultralight scalars it still make sense to discuss their effect as a coherent gravitational 
waveQ 

In the opposite limit of large mass, the single quanta of the scalar field will behave as 
particles rather than waves and will interact incoherently with the detector, just with separate 
hits and since they only interact gravitationally they will be totally unobservable. 

Consider the action describing the Jordan-Brans-Dicke theory (|]), with in addition a 
potential term for the scalar: 



Sp + Sr< 



Smbfcm, 9nv] = /d 4 aV = 5£ m [V>m,ff At i/] • (59) 



In analogy with the massless case, to obtain the field equations we vary the action S with 
respect to the metric g^ v and to the scalar ip\ then, we linearize the equations near the 
background (rj^jipo), where (po is a minimum of V. The equations of motion in vacuum 
(where the matter energy-momentum tensor is zero) are 

n^-^n^ = -d^Q + n^m (eo) 

□ <J>-m 2 $ = 0, (61) 

where 



3 Of course, in general such a ultralight scalar would mediate an unacceptably strong fifth force: Newton's 
law is tested down to length of order ~ cm, so that an hypothetical fifth force could be mediated only by 
a particle with mass m > (1 cm) -1 ~ 2 x 10 -5 eV (see e.g. ref. Q). Here, however, we are considering 
Brans-Dicke theory, with ljbd > 600: the corrections to Newton's gravity, in the weak field limit, are smaller 
than 0(1/wbd), so that even a massless scalar is compatible with experiments. 
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1Z^ V ,1Z are the linearized (with respect to h^ u = (g^ — rj^) ) Ricci tensor and scalar 
curvature; 

$ = -(if - ip )/lfo ; 

2 _ V"((f ) ip 



m 



(3 + 2w B d) 



Again in analogy to the massless case, the linearized Riemann tensor IZ^vpa and the equations 
(|60| ) and (|6l]) are invariant under gauge-transformations 

f VW = (62) 

I ¥>(a0 -» p'(x) = <p(x) ; 

we then define 

V - ^ f J/» - (63) 



and choose the Lorentz gauge (in analogy to electromagnetism) 

V" = , (64) 

by means of a transformation such that Ue v = d^9^ u . In such a gauge the field-equations 
become wave-equations: 

□V = ( 65 ) 

= m 2 $ , (66) 

whose solutions are plane waves (and superpositions of them): 

M*) = VW^*" + c.c. (67) 

= b(k)e iqaxa + c.c. , (68) 



where: 



(wo;*;) w = A:ee|A;| _ 

i 2 ■ 



(a> m ;A:) w m = 



as one sees from fl69[), the dispersion law for the modes of $ is that of a massive field, and 
the group- velocity of a wave-packet of $ centered in A; is 

v = -zr = — ' ( 70 ) 

OK w m 

exactly the velocity of a massive particle with momentum k. 

As in the massless case, we have a residual gauge-freedom: we remain in the gauge fl64{ ) 
by transformations with De u = 0. Notice that the Lorentz gauge is a transversality condition 
on the field 9^ u 

= => k^ lu = (71) 
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but it says nothing about the transversality of the perturbation h^ u = 9^ u — rj^ ( ^9 — $ 



And, in fact, it is easy to show that, due to the field equation (66), it is impossible to choose 
a gauge analogous to (g). In fact, as in the massless case discussed in app. A, to choose such 
a gauge, we should operate a transformation with: 

J Ue^ = 

1 = -Jfl + fc, (72) 



but now this is not compatible with the wave equations (65), (|36|). Analogously, we cannot 
impose any linear relationship between 9^ and $: so, for example, we cannot choose a gauge 
in which hnQ = 0, which in the massless case instead permits to integrate immediately the 
geodesic equation. 

Instead, we can choose a gauge in which a gravitational wave propagating in the +z 
direction takes the form: 

M*, *) = A{+ \ 1 ~ z ) ffi + ^\t - z) e W + z) Vfll/ , (73) 

where effi and e]L* are the polarization tensors @. This is similar to the gauge choice 
discussed in app. B, but now, because of the non-linear dispersion law, $ depends separately 
on t and z, instead of depending only on their difference t — z. This makes the discussion of 
the physical effect of the wave different from the massless case. 

To study the effect of the scalar wave on test-masses, we now make use of the proper 
reference frame of one of these masses (as the beam-splitter of an interferometer, for example): 
the analysis thus can be performed in newtonian terms (see sect. 13.6]), in the sense that 
the spatial coordinates x 3 represent proper distances for the observer 'sitting' on the beam- 
splitter, the time variable t represents his proper time, and the effect of the gravitational 
wave on test-masses is described by the equation of motion: 

-JL = -& oko x k , (74) 

where the quantities R J oko are the so-called electric components of the Riemann tensor. The 
problem is thus reduced to calculate i? J ofco m the proper reference frame of this observer: 
this is connected by an infinitesimal coordinate transformation to the reference frame where 
eq. (|73|) holds, and since the linearized Riemann tensor TZ^ pa is invariant under infinitesimal 
gauge transformations, we can compute it from the metric eq. (|73|). Restricting to a purely 
scalar GW, a straightforward computation gives 

-d 2 t 

ft*c«o = \ | -dj | $(t,z) = —^Tij ~L id m*$ , (75) 

m 2 

where = (5ij — hihj) is the transverse projector with respect to the direction of propagation 
n, and Lij = hihj the longitudinal projector. The equation of motion (|74]) thus becomes: 



d 



2 



dt 2 x i [ k d t® ) ■ T H x j + ( - ) ■ Lij x j , (76) 
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from which one immediately sees what is the effect of the dilaton mass: it generates a 
longitudinal force (in addition to the transverse one) proportional to m 2 . In the limit m — > 
one recovers the treatment of the massless case, confirming again, also from the point of view 
of the proper reference frame, the results of sect. 2. 

To better understand the contribution of m, it is convenient to restrict to a wave packet 
<3?(i, z) centered at a frequency uj m = \Jk 2 + m 2 ; in this case 

= (-lufa&j • (t^ + ^- Li^j Xj . (77) 

From this equation the physical effect of the mass can be read quite clearly. In particular 

- for ultra-relativistic momenta (u m ~ k ^> m) the longitudinal components of the force 
becomes negligible with respect to the transverse one; 

- in the non-relativistic limit tu m — * m and the stress induced by the wave becomes 
isotropic, as Ly + = Sy. 

Now one can compute the pattern function of a detector for these massive scalar waves. 
In a generic gauge the response of the detector is obtained contracting the detector tensor 
Dij with the tensor ffojo, which in the massive case is proportional to co^Tij - s r r n 2 Lij. Note 
that the pattern function of an interferometer for detection of such massive scalar waves is 
identical to the one computed in the massless case: in fact, the contraction of w^Tjj + m?Lij 
with is the sum of a term containing DijTij and a term containing DijLij; both these 
terms are proportional to hihj, since is traceless, for an interferometer, and, so the 
only dependence on the mass m is in the overall factor of the signal and is independent of 
n. In particular, for tu m — > m the signal goes to zero. The pattern function is instead mass- 
dependent in the case of a detector with a non-traceless response tensor, as for the monopole 
mode of the sphere, or for the common mode of an interferometer, for which = UiUj+ViVj. 



5 The string dilaton and moduli fields 

It is clearly important to understand whether these results, that have been obtained in the 
context of Brans-Dicke theory, can be applied to the physically more interesting case of the 
dilaton and the other scalar fields predicted by string theory. 

While the dilaton-graviton sector of the low energy sector of string theory is the same as 
a Brans-Dicke theory, the situation is quite different for the interaction of the dilaton with 
matter. Since in the string case wbd = — 1, the dilaton is coupled with a strength of the same 
order as the graviton, and produces unacceptable deviations from general relativity, unless a 
non-zero dilaton mass is generated see e.g. pj. The radius of the non-universal force that it 
mediates must be smaller than about 1 cm, or m > 2 x 10 -5 eV. Therefore the analysis of the 
previous section, which was valid for m < 4 x 10 -12 eV, does not apply to a massive dilaton. 
Indeed, it is also in general not easy to reconcile such light scalar particles with cosmology, 



see e.g. [21, 22], although there are mechanisms that solve the cosmological problems created 
by light scalars, typically introducing a second short stage of inflation that dilutes the dilaton 
overproduced by oscillations around their quadratic potential. 
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Actually, there is the possibility to circumvent this bound on the mass, with a mechanism 
that has been proposed by Damour and Polyakov ||]. Assuming some form of universality 
in the string loop corrections, it is possible to stabilize a massless dilaton during the cos- 
mological evolution, at a value where it is essentially decoupled from the matter sector. In 
this case, however, the dilaton becomes decoupled also from the detector, since the dimen- 
sionless coupling of the dilaton to matter (a in the notation of ||) is smaller than 10~ 7 (see 



also [23]). Such a dilaton would then be unobservable at VIRGO, although it could still 
produce a number of small deviations from General Relativity which might in principle be 
observable improving by several orders of magnitude the experimental tests of the equivalence 
principle [24|. 

So, in both cases, the analysis done for Brans-Dicke theory does not appear to be relevant 
for string theory. However, it is clear that our present understanding of the string dilaton 
and moduli is incomplete and presents a number of unsettled issues, including the non- 
perturbative mechanism for mass generation, or the stabilization at the minimum of the 
potential [25], and a definite conclusion is probably premature. 
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A Gravitational waves in the Jordan— Brans— Dicke theory 

The Jordan-Brans-Dicke theory is described in the Jordan-Fierz frame by the action 

S = S g [g^, if] + S m [V> m , diiv\ (78) 

(79) 



1 ',4„ 



s m = [d A %V-g £mbl>m,9nv] ; (so) 



where: 



• guv is the metric tensor, with which one constructs all the covariant quantities, such as 
the scalar curvature R, covariant derivatives V, etc.; 

• tp is a scalar field; 

• wbd is a parameter; 

• ip m are the matter fields, such as fermions, gauge fields, etc.; 

• C m is their lagrangian. 

To obtain the field equation one as to vary the action S with respect to g^ v and to (p; this 
leads to 



d^ip dyip - \g^g af5 d a <p dpcp 



+ 
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+- 



r 



V^^f/j - g^g^ Vadpip 



r 



where 



T = 



2 5S m 

r^jSg^ 



(81) 
(82) 

(83) 



is the matter-fields energy- momentum tensor. 
Eq. (gl]), multiplied by g^ v , gives 



that, substituted in eq. (82), leads to: 

(3 + 2lo bd ) 



(84) 



(85) 



Eqs. (0) and (||) are our basic field equations. We now study small perturbations around 
a background configuration: we choose as background the Minkowski metric and (p{x) = 

(fQ = - — — — ——r^ , in order to have the correct post-Newtonian limit and to obtain General 
(3 + 2wbd)& 

Relativity when lobd — > oo , as showed in pq] ; we consider 



<p(x) 



rj^u + h pu (x) 
ip + 5ip{x) , 



(86) 



with \h pu \ <C 1 and \5ip\ <C (fo- We call lZ^ upa , lZ pu and 1Z the linearization to first order in 
hfju, of the corresponding quantities R pU pa-, Rp,v an d R] one has |2C] 



Ti^upa = f {dpdpKv + dvdchrf - d a dph pu - d^Kp} . 



The linearization of the field-equations (^T]) and (pq) in vacuum (T pu = 0) gives 



K 



pu 



0, 



(87) 

(88) 
(89) 



with <I>(x) = —5ip(x)/ipQ. In analogy to General Relativity, we can define a transformation 
acting on our fields that leave unchanged the linearized Riemann tensor lZ^ upa (and conse- 
quently equations (]88|) and (p9|)): such an infinitesimal gauge transformation with parameter 



e>" is 



(90) 



It is straightforward to verify, by direct substitution in (|87|), that lZ^ upa is gauge-invariant. 
We want to use this gauge freedom to obtain a wave-equation. We thus define 



' pu 



h 



pV 



= rre^ 



-h + 4$ , 



(91) 
(92) 
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with h = rf ,v h IJa /', the transformation that expresses hn U in terms of 0u U has the same form 



/V = ^ ~ ^Vtiud + r)^ v § (93) 
h = -6 + 4$ . (94) 

Substituting ( |93"| ) into ( |38| ) one obtains the field equation for 

□V - (d a e av ) - d u (d a e a ^) + (d a e aP ) = o . (95) 

Inserting eq. (|90|) into eq. (EJT1) one finds immediately that under gauge-transformations 



0(z) 
<3>(x) 



6{x) + 25° e a 



By choosing such that Oe u = d^Onvi one nas 



and, so, a wave equation for 0' (see eq. fl95|)) 



(96) 



(97) 



(98) 



In this gauge (we will call it Lorentz gauge, in analogy to electromagnetism) the solutions are 
plane waves and their superpositions (we omit the ( )') 



<£(x) 



V(*0 e ^ + c - c - (") 

b(k)e ikax <* + c.c. , (100) 

with the following conditions (deriving from the field-equations and from (j97|)) 

(101) 

, /( , ; (102) 



the latter is a transversality condition for A^ u . 

Once we have chosen the Lorentz gauge, we can still operate transformations with De^ = 0; 
we thus take e M such that 

f De" = 



(it is possible bacause in vacuum in Lorentz gauge □# = □<& = 0); 
one has (see ([)(]) and (|93|)) 

= 2$ => v = V ; 

that means that /t„ v too is a plane transverse wave. 

Again, we have not yet completely fixed the gauge: we satisfy our conditions 



(103) 



(104) 





2$ 



(105) 
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also by operating gauge transformations with 

f a e M = n 

{^ = 0. < 106 > 

Consider the case in which the wave is propagating in +z direction: then 

k^ = (A;, 0,0, k) (107) 

fc^V = A 0u = -A 3v (108) 

A u0 = -A v3 (109) 

A 00 = ~A m = +^33 • (HO) 

Let us make a degrees-of- freedom counting for A^ v : 

- we started with 10 (= independent components of a symmetric tensor); 

- transversality leads to 7: it "kills" only 3 (instead of 4) because of symmetry of A^; 

- the condition 8 = 2<3? leads to 6: we choose Aqq, An, A22, A21, A31, and A32 as 
independent components; 

- further gauge freedom permits us to put to zero 3 of the 6 components (three rather 
than four, because of condition = 0). 

Thus taking 

f e„(x) = e^k) e ikax <* + c.c. . * 

1 = K 1 



the action of the gauge-transformation on A^ u is (see (|9q) and (|99f)) 

Ap, v — > = .A^i/ — ik^e^ ; (112) 
or, for the 6 components we are interested in, 



(113) 



Notice that An, A22 and A21 = A\2 are invariant: we thus choose eo, ei, £2 in order to "kill" 
the others. We have now completely fixed the gauge. 

Does h,j,u depend on the field $? Eq. ( |l04l) tells us h = h\\ + ^22 = 2$. 
Summarizing, we can say that, zn i/iis gauge, the metric perturbation /tuj, produced by a 
gravitational wave propagating in the +z direction takes the form 

V(* - *) = ^ {+) (* " z) e<+> + A( x )(i - z) eg + *(* - z) eg , (114) 



A 00 - 


-> Aw + 2 £&e 


An - 


- An 


A22 - 


■* A 22 


A21 - 


■+ A 2 i 


A31 - 


-> A31 - ifcej 


A 32 - 


-» A 32 - i/ce 2 . 



where e+^ x ' s are given in eqs. (0) and 
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B Equivalence of the two gauges 



In this appendix we show that it is possible to choose a gauge (which we denote writing a 
tilde over quantities evaluated in this gauge) in which the metric perturbation has the 
form (IIT41), with AM 



AM, AM 

;(s ' ) = diag 



A(*\$ = $,eM 



J+) zM 



,M 



, but 



-1,+1,+1,+1) =r)»u ■ 



(115) 



The quantities without a ( ) refer to the gauge used in the previous appendix and in section 
2.1. We call "transverse" the former gauge, and "conformal" the latter: we want to find a 
gauge transformation that passes from the transverse gauge to the conformal one. In order 
to have a wave equation for 6^ one has to keep the Lorentz gauge condition: this condition 
is mantained by imposing rJe» = 0. By choosing 







(116) 



which is possible because in Lorentz gauge in vacuum n6 = 0, we obtain a traceless 6' and 
so (see eq. fl93|)) 

Ku = 0' + $ . (117) 



In analogy with appendix A, acting with a transformations with 

1 Ue» 








(118) 



we can now eliminate the appropriate components of 9' (see equations flll3l )), in order to 
obtain exactly L v . It is then easy to prove that AM = AM an d AM = AM j by examining 
the action of gauge transformations on those amplitudes. 

The two gauges are thus equivalent: and in fact it is easy to exibit the coordinate trans- 
formation that relates them (the coordinates with a ( )' denote the transverse gauge and we 
now limit to the case of a purely scalar wave) 



\l{t - z) 
t+ l -I{t-z) , 



;ii9) 



where 



In fact 



I{t-z) 



t-z 



(120) 



' dx' = 


dx 


dy' = 


dy 


< dz' = 




dt' = 


I 



(121) 



<f> dz 
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and 



ds 2 



(1 + <f>){-dt 2 + dx 2 + dy 2 + dz 2 ) = 77^ dx m dx' v + 3> (dx' 2 + dy' z ) 



./2n 



(122) 



As we have seen in section 2.1, the physical meaning of the primed coordinates is that they are 
comoving with free-falling test-masses, initially at rest, and t' is proper time. We can check 
this assertion in the conformal gauge, by making use of the solution of the geodesic equations 

By substituting eq. (22) into the transformations (|Tl1 



of motion found in |15| , see eq. (| 
one finds 



At) 
y'it) 
At) 
fit) 



x 

JJO 

t + z{t) 



(123) 



zq = r(t) + const 



that proves that, in the primed coordinates, bodies initially at rest remain at rest. 



C Relationship between Sh(f) and Qip(f) 

When dealing with stochastic backgrounds of ordinary GWs one defines Sh(f) as follows [18]. 
One expands the metric perturbation in plane waves 

/+oo r 
df / dtt fi h A (f, n) e 2m eg.(fi) , (124) 

where the ensemble average of the Fourier modes is 

( h* A (f,n) h AI (f',n') )=5(f- /') i- ,5 2 (n,n') 5 AA , \ S h (f) , (125) 



and <5 2 (n, n') = 5(cf> — (f) 1 ) 5(cos9 — cos9') . By inserting eq. (124) in the expression p gw = 
l/(327rG) ( hijhij ) for the energy density of the background, one obtains 

Pgw = ^~ + °° d(1 ° g /} f (27F/)2 Shif) ' (126) 

so that 

W) = ^js %(f) ■ ( 127 ) 
In the case of scalar waves, we write 

4ht,x) = [ + °°df fdn & h {B) (f,n) e *rf e g.)(fi) , (128) 

J -co J 

and 

( h* {s) (f,h) h {s) (f',h>) ) = S(f- f) ±- J 2 (n,n') \ S h (f) . (129) 

We want to relate this new function Sh(f) to the energy density p v of the field ip. It is 
convenient to rewrite the field-equation of the Jordan-Brans-Dicke theory in the Einstein 
frame, that it is related to the Jordan-Fierz one by the conformal transformation 

9% = (^) 9^ • (130) 
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In this frame, the field-equation for the metric has the standard General Relativity form 



„E _ I E nE ^ T E , 3 + 2(J BD 



d^p d u (p - \g^ y 5e d a ip dpip 



and the energy-momentum conservation law becomes 
lE (3 + 2u; B d)<A) 



(131) 



VS 



- /if 



167T (/? 2 



(9 /t (/J <%l/2 - fl^f <9 a (/7 9/3^ 



vs. 



T E -I- T^ 3 ) 

/if ' ■*■ flU 



(132) 



so we define T^/ as the energy- momentum tensor of the field </3. To first order in h E 



{9% - Vtiv) and $ 



1 - <fo 



we have 



(133) 



the energy density is 



= T (rf _ ( 3 + ^BD^O 

/V — J oo — ' 



32vr 



$ 2 + 



(v*) 2 ' 



(134) 



that, by using field equation = and averaging over several wave-lengths, reduces to 

(4 + 2a; BD ) , ^ 2 



32vrG 



-(2^ ) 



(135) 



where we used (po = (4 + 2u;bd)/G ! (3 + 2ubd), as discussed in app. A. As h ( ff = $eg J , using 



eqs. ( |128| ) and ( |129|) we have 



32vrG « « 
2(4 + 2l^bd) f f=+0 ° 



so that 



327rG "'/=o 



s h (f)- 3H2 ° 



d(logf)f(2nf) 2 S h (f) 



1 



4vr 2 /3 2 + wbd 



(136) 
(137) 

(138) 
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